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■ ABSTRACT 

> ■ 

Conservation laws of the nonlinear Schrodinger equation are studied in the presence of 
higher-order nonlinear optical effects including the third-order dispersion and the self- 
steepening. In a context of group theory, we derive a general expression for infinitely 
many conserved currents and charges of the coupled higher-order nonlinear Schrodinger 
equation. The first few currents and charges are also presented explicitly. Due to the 
higher-order effects, conservation laws of the nonlinear Schrodinger equation are violated 
in general. The differences between the types of the conserved currents for the Hirota 
and the Sasa-Satsuma equations imply that the higher-order terms determine the inher- 
ent types of conserved quantities for each integrable cases of the higher-order nonlinear 
Schrodinger equation. 
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In the ultrafast optical signal system, the higher-order nonlinear effects such as the 
third-order dispersion, the self-steepening, and the self-frequency shift become important 
if the pulses are shorter than T < 100/s [[[J. The use of optical pulses with distinct polar- 
izations and / or frequencies also require the consideration of nonlinear cross-couplings be- 
tween different modes of pulses. Inclusion of both the higher-order and the cross-coupling 
effects lead to the study on the coupled higher-order nonlinear Schrodinger equations 
(CHONSE) which are not in general integrable except for special cases of coupling con- 
stants. Those integrable cases of coupling constants have been classified in association 
with Hermitian symmetric spaces |J. It is also well known that soliton equations which 
can be integrated by the inverse scattering method possess infinite number of conserved 
quantities. For example, the nonlinear Schrodinger equation (NSE) has infinite number of 
conserved charges in addition to the ones corresponding to the energy and the intensity- 
weighted mean frequency. However, the effect of the higher-order and the cross coupling 
terms on the conservation laws has not been considered up to now. 

In this paper, we make a systematic study on the conservation laws in the presence 
of the higher-order and the cross-coupling terms. We first indicate that except for the 
energy conservation, other conservation laws of the NSE such as the conservation of the 
intensity-weighted mean frequency do not hold due to the higher-order effects any more, 
unless the higher-order terms are of a unique type. In the case of integrable CHONSE, we 
derive general expressions of infinite number of conserved currents and charges from the 
Lax pair formulation utilizing the properties of the Hermitian symmetric space. From the 
general expressions, explicit forms of the first few conserved currents and the associated 
charges of the Hirota and the Sasa-Satsuma equations are calculated. We then explain 
the correlations of conservation laws between the two integrable cases of the higher-order 
extension of the NSE. 

In order to illustrate the issue, we first consider the NSE including the higher-order 
terms. In a mono-mode optical fiber, the propagation of a ultrashort pulse is governed 
by the higher-order nonlinear Schrodinger equation || 

dip = i( 7 i0V + + 73<9 3 V> + 74WrV) + 75<9(H 2 )^, (1) 

where 8 = d/dz and d = d/dz are derivatives in retarded time coordinates (z = x, z = t — 
x/v), and if) is the slowly varying envelope function. The real coefficients 7$ (i = 1, 2, 3, 4) 
in the first four terms on the right hand side of Eq. (TO) specify in sequence the effects 
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of the group velocity dispersion, the self-phase modulation, the third order dispersion, 
and the self-steepening. With appropriate scalings of space, time, and field variables, one 
can readily normalize Eq. (|J) so that 71 = 1, 72 = 2, 73 = 1 which we assume from 
now on. The remaining coefficient 75 in the last term is complex in general. The real 
and the imaginary parts of 75 are due to the effect of the frequency- dependent radius 
of fiber mode and the effect of the self-frequency shift by stimulated Raman scattering, 
respectively. It is well known that the above equation becomes integrable if 74 = —75 = 6 
(Hirota case) J| or 74 = —275 = 6 (Sasa-Satsuma case) f|. In the absense of higher-order 
terms (73 = 74 = 75 = 0), Eq. ([!]) possesses infinite number of conserved charges among 
which the first three charges are 

/CO 
\tfj\ 2 dt , 
-00 

/oo 
-00 

/oo 
(d^diJ - \ij\ 4 )dt, (2) 
-00 

where Qi represents conserved energy, and Q2 the mean frequency weighted by the inten- 
sity of optical pulses. In the conventional NSE where the time and the space coordinates 
are interchanged, Qi,Q2 and Q3 respectively correspond to conserved mass, momentum 
and energy. If we include higher-order terms, Qi are not necessarily conserved but subject 
to the relations; 

BQ 1 = 0, 

/oo 
c#| 2 , frW - Wtydt 
-00 

/oo 
d\il)\ 2 d^*d^dt . (3) 
-00 

The calculations indicate that the charge Q\ which corresponds to energy is conserved for 
all values of 74, 75 while Q2 and Q3 are conserved provided 74 + 75 = and 374 + 275 = 6, 
respectively. Note that Qi and Q3 are conserved simultaneously only for the specific value 
74 = —75 = 6 that is precisely the Hirota case. It is interesting to observe that integrability 
does not always imply the same types of conserved currents in the presence of higher-order 
terms. Another integrable case of the Sasa-Satsuma equation, where 74 = —275 = 6, in 
fact does not have Q 2 and Q 3 in Eq. (|3|) as the conserved charges. This consequence is 
rather remarkable in view of the fact that integrable equations possess infinite number of 
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conserved quantities. We will show, however, the Sasa-Satsuma equation also possesses 
infinitely many conserved charges of different types other than the ones of the Hirota 
equation. 

In case we include both the higher-order and the cross-coupling nonlinear effects, the 
propagating system is governed by a CHONSE. Without understanding physical settings, 
it would be meaningless to write down any general expression of the CHONSE. However, 
as explicitly derived in 0, there exists a group theoretic specification which admits a 
systematic classification of integrable cases of the CHONSE. In the following, we consider 
a group theoretic generalization of the NSE and define the CHONSE in association with 
a Hermitian symmetric space. By solving the linear Lax equations iteratively, we derive 
infinite number of conserved currents and charges for the CHONSE. For the later use, now 
we briefly review the definition of Hermitian symmetric spaces f?|, H and the generalization 
of the NSE 0, |9| according to the Hermitian symmetric spaces. 

A symmetric space is a coset space G/K for Lie groups G D K whose associated Lie 
algebras g and k, with the decomposition: g = k © m, satisfy the commutation relations; 



A Hermitian symmetric space is the symmetric space G/K equipped with a complex 
structure. One can always find an element T in the Cartan subalgebra of g whose adjoint 
action defines a complex structure and also the subalgebra k as a kernel, i.e., k = {V G 
g : [T, V] = 0}. That is, the adjoint action J = adT = [T, *] is a linear map J : m — > m 
that satisfies the complex structure condition, J 2 = —I, or [T, [T, M]] = —M for M e m. 
Then, we define a CHONSE as| 



where E and E = [T, E] are extended field variables belonging to m. The arbitrary 

constant a may be normalized to 1 by an appropriate scaling but we keep it in order 

to exemplify the higher-order effects. Also the cross-coupling effects between different 

modes of polarizations or frequencies are accommodated in the matrix form of E which is 

determined by each Hermitian symmetric space. For example, in the case where G/K = 

4 We restrict to symmetric spaces AIII, CI and Dili only so that the expression of 
CHONSE becomes simplified @. 



[k, k] C k, [m, m] C k, [k, m] C m 



(4) 



dE = d 2 E - 2E 2 E + a(d 3 E + (3 x E 2 dE + (3 2 dEE 2 ) 



(5) 
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SU(N + l)/U(N), matrices E and T are represented as 



E 



( o Vi 
o 



V— ^ 



TV 









o J 



T 



\ 




-\) 



(6) 



and the CHONSE becomes the higher-order vector nonlinear Schrodinger equation, 



N 



a 2 ^ + 2(EI^I 2 )^ 



i=i 



— a 



AT TV 

/ME \H 2 Wk + /3 2 (E - ^ 

i=i i=i 



fc = 1,2, 



N. 



(7) 



This equation is an obvious generalization of Eq. (|l|) to the multi-component case. It is 
easy to see that Eq. (|7]) with N — 1 and Pi = /3 2 = —3 is precisely the Hirota equation. As 
verified in 0, Eq. (|) is integrable if Pi = p 2 = — 3 because in such a case the CHONSE 
admits a Lax pair. That is, Eq. (j5|) with Pi = P2 — — 3 arises from the compatibility 
condition ( L g ] = 0) of the associated linear equations, 

L z fy = [d + E + AT]* = 0, 

L g tf = [B + U° K + U° M + \(U 1 K + U 1 M ) + \ 2 (U 2 M + T)-a\ 3 T}^ = 0, (8) 

which holds for all values of the spectral parameter A. The entities U l K and U l M in L 2 are 
given by 



U° K 



-EE - a[E, dE], U° M = BE + a(d 2 E - 2£ 3 ), 
aE£, = E — adE, U 2 M = -aE. 



(9) 



Here, the subscripts K and M signify that they belong to the subalgebra k and the 
remaining complement m, respectively. The algebraic decomposition can be also extended 
to a more general case including the matrix solution * = * % + \&m with the properties 
that [T, ^k] = 0, [T, £ m, and the following multiplication properties; 



[T, tf]^] = [T, *i,*3,] = 0, [T, ^ 2 M ] G m 



(10) 



The adjoint action of the element T in the Cartan subalgebra together with the complex 
structure condition, if applied to the decomposition, lead to a couple of general identities 
for any M 1 , M 2 G m: 



[T, M X M 2 \ = M\M> + MiM 2 = 0, M;M> = M X M 2 



These identities are useful for many calculations, for example, in deriving conserved cur- 
rents or in verifying that the CHONSE in Eq. (|5]) is equivalent to the compatibility 
condition of the Lax pair in Eq. (|||). 

Having presented necessary ingredients, we are now ready to derive infinitely many 
conserved currents and charges of the integrable CHONSE by solving the associated linear 
equations in Eq. (|8|). In order to make use of the algebraic properties of Hermitian 
symmetric spaces, we make a change of the variable \1/ in Eq. (|5|) by 

$ = ^ exp[(A^ + (A 2 - aA 3 )^)T], (12) 

which results in the change of the multiplicative term T^f to the commutative term [T, <3>] 
in the linear equations. The adjoint action, [T, <&], allows the splitting of the linear 
equations for $ into the K- and the M-components as explained below. Let us first 
assume that the linear equations can be solved iteratively in terms of 

OO 1 

*{z,z,\) = £ T^(*k(M) + SJtrCM)), (13) 

n=0 A 

where and $^ denote the decomposition of a coefficient $ n satisfying the properties 
in Eq. (0). Then, the n-th order equation (n > 0) separates into the K- and the 
M-components such as 

d<& n K + E$ n M = 0, (14) 
dn i + E$ n K + [T,$ n A f 1 ] = 0, (15) 

while the <9-part of the linear equation becomes 

mi + u° K ^ K + u° M ^ M + u 1 ^ 1 + umt 1 + umr = o, (i6) 

M n M + U° K n! + 02r*x + Uhnt + Ul^ n K +1 + Ul^T + [T, - a[T, nf] = 0- 

(17) 

In addition, there are equations arising from the positive powers of A, which can be given 
by Eqs. provided that n = -1, -2, -3 and $^ <0 = = are defined. These 

equations can be solved recursively for (n > 0) starting from a consistent set of 

initial conditions; 

$° M = 0, & K = -iI, & u = -iE. (18) 
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Note that Eq. flTSP can be solved for by using the complex structure condition. That 
is, = — [T, [T, = [T,<9$^] + Thus, ^Ij 1 is obtained directly provided 

that $^ and are determined. Contrary to which is obtained algebraically, 
can be obtained by a direct integration of Eq. ([14]). In fact, is over-determined due 
to Eq. ( |T6| ) as well. Thus, in order for to be integrable, the compatibility condition 
that [<9, 5]$^- = should be required. The condition is satisfied provided the integrable 
CHONSE holds. In this case the compatibility condition gives rise to infinitely many 
conserved currents labeled by integer n such that dJ 1 ^ + dJ 1 ^ = 0; 

Jl = -d$ n K = E<S> n M (19) 
Jl = d$ n K = -(dE + ad 2 E - 3aE 3 )$ n M - aEd 2 $ n M - (E - ocdE)[T, d§ n M ] (20) 

In order to derive the local currents explicitly, we solve the recurrence relations in Eqs. 
flnj)-(O) with the initial conditions as in Eq. (13). The first few conserved currents are 
listed below; 

J l K = -iEE, 

J\ = -i[E,dE]+ta([d 2 E,E] + dEdE -3E 3 E), (21) 

for n — 1, and 

= -idQ^&x + iEdE, 
= +id® 1 K <f> 1 K - i(Ed 2 E + dEdE - E 3 E) 

-ia(Ed 3 E + [d 2 E, dE] + dEE 3 - 2E0EE 2 - E 2 dEE - AE 3 dE), (22) 
= -i{d<& 2 K <& x K + d® 1 ^ - id&K&x&x) + i[Ed 2 E - E 3 E), 
= i(d<S> l K <S> 2 K + d$ 2 K <$> l K - id&K&K&x) 

+ i(Ed 3 E - dEd 2 E + dEE 3 - 2E8EE 2 - E 2 dEE - 2E 3 dE) 
- ia(Ed 4 E + d 2 Ed 2 E + dEd 3 E - 5E 3 d 2 E - E 2 d 2 EE - 3Ed 2 EE 2 
+d 2 EE 3 - 2dEdEE 2 - dEEdEE - 2dEE 2 dE - 3E8E8EE 
—5EdEEdE — 3E 2 dEdE + AE 5 E), (23) 

for n = 2 and n = 3, respectively. Note that currents J% and J\ for n > 2 contain 
non-local terms $^ with m < n. Fortunately, these non-local terms can be separated 



J' 2 



K 



T 2 



T 3 



T 3 
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from the conservation law if we consider a scalar expression of the conserved current by 
taking an appropriate trace as follows; 



S n K = Tr(PJZ), Sl = Ti{PJ n K ) (24) 

The parameter P is any matrix entity which commutes with matrices or we may 
choose P = C\l + c{P for arbitrary constants C\ and c<i- For instance, we have for n = 2, 3 

S 2 K = -d(TrP{ l -(^ K ) 2 }) + TrP(zEdE), 

S 2 K = <9(Tr P{ l -{^ K f}) 

+ Tr P{-i(Ed 2 E + dEdE - E 3 E) - ia(Ed 3 E + [d 2 E, BE] - 6E 3 dE)}, (25) 
S 3 K = -9(TrP{^($^-i($i0 3 )})+TrP{z(^ 2 E-^ 

S 3 K = a(TrP0($^-i($^) 3 )}), 

+ Tr P{i(Ed 3 E - dEd 2 E - AE 3 dE) - ia{Ed 4 E + d 2 Ed 2 E + dEd 3 E - 8E 3 d 2 E 
+2d 2 EE 3 + E 2 dEdE - dEEdEE + dEE 2 dE - 5EdEEdE + AE 5 E)}. (26) 

The derivations show that the non-local terms appear as total derivative terms thus they 
are conserved separately. Dropping the non-local terms and integrating over the time 
coordinate, we obtain infinite number of global charges which are conserved in space, i.e. 
dQ\ = where 

/+oo 
dtS n K . (27) 
-oo 

For the case of G/K = SU (N + 1)/U(N) as mentioned in Eq. (|B|), we work out explicitly 
and obtain the conserved charges 

+oo 



N 

Qk= I rft E^fc, (2? 
fc=i 



for n = 1 and 



r +oo N 

Q 2 K = i dt^tdA - Wi), (29) 
J -°° k=i 

Qk = / dt[J2 a^w* - (E rM% (so) 
J -°° k=i k=i 

for n = 2 and n = 3, respectively. Conserved charges for other cases of integrable 
CHONSE can be similarly obtained from the specification of E and T as classified in H. 



As noted in Eq. (|3|), the types of charges Q2 and Q3 are not conserved in the Sasa- 
Satsuma case where 74 = —275 = 6. Nevertheless, the Sasa-Satsuma equation equivalently 
possesses infinitely many conserved charges of different types as well. This seemingly 
contradicting characteristics can be explained by the fact that the Sasa-Satsuma equation 
arises from the discrete ^-reduction of the SU{2>)/U{2) CHONSE combined with a point 
transformation 101. In this case, matrices E and T can be denoted as 




E = \ -ib* , T 



(\ 



0-1 |. (31) 

Vo -\ 



Since the charge Q\ in Eq. (|27|) is invariant under the point transformation, we can 
also calculate the first few conserved charges of the Sasa-Satsuma equation using the 
expressions E and E = [T, E] given in Eq. fl3~l|) . The resulting charges of the Sasa- 
Satsuma equation are 



q) { = / dtr^ 
j —00 

Qk = 0, 

r+00 

Q\ = / dt[3dib*dib -6(ip*ij) 2 -i(ib*dik - dib*ik)}. (32) 



If the charges in Eq. (^) are compared with those of the Hirota type in Eq. (§) (or 
equivalently Eq. (p8|)-(p0|) for N = 1), we note that the current for n — 1, which corre- 
sponds to energy, is the same but other currents are of different types. Remarkably, in 
Eq. (^) the current for n = 2 turns out to be trivial while the current for n = 3 is a 
new type that is seemingly combination of currents for n = 2, 3 in Eq. (^). From Eq. 
([!]) with normalized coefficients 71 = 72/2 = 73 = 1, one can readily confirm that the 
current = 3dip*dtp — 6(i/j*ip) 2 — i(ib*dib — dip*ip) is conserved only if 74 + 75 = 3 and 
374 + 275 = 12. Solving the equations results in 74 = —275 = 6 that definitely leads to 
the Sasa-Satsuma case, to be compared with Eq. (^) for the Hirota case. 

To summarize, using the properties of Hermitian symmetric space we have constructed 
the Lax pair formalism of the coupled higher-order nonlinear Schrodinger equation and 
derived a general expression of infinite number of conservation laws. Remarkably, the 
conserved currents and charges for both the Hirota and the Sasa-Satsuma equations 
are calculated from the general expressions, accompanying the reduction procedure. We 
have shown that, except for the Hirota case, the current conservations of the nonlinear 
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Schrodinger equation are in general broken by the higher-order effects. The types of con- 
served currents and charges for the Sasa-Satsuma case are different from the types for 
the Hirota case except for the energy conserved irrespective of all the higher-order effects. 
These differences may leave scope for more physical explanations and applications in the 
further study of higher-order effects including numerical analysis. 
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